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Waveguides for walking droplets 
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When gently placing a droplet onto a vertically vibrated bath, a drop can bounce without coa¬ 
lescing. Upon increasing the forcing acceleration, the droplet is propelled by the wave it generates 
and becomes a walker with a well dehned speed. We investigate the conhnement of a walker in dif¬ 
ferent rectangular cavities, used as waveguides for the Faraday waves emitted by successive droplet 
bounces. By studying the walker velocities, we discover that Id conhnement is optimal for narrow 
channels of width of D ~ 1.5Af. We also propose an analogy with waveguide models based on the 
observation of the Faraday instability within the channels. 
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INTRODUCTION 


When a droplet is gently placed at the surface of a vibrated liquid bath, inhibition of coalescence can be obtained 
depending on the forcing parameters [1]. Thanks to the air layer between the drop and the bath which apply a 
lubrication force, the contact can be avoided. The subsequent dynamics is rich: numerous bouncing modes [2], anti¬ 
resonance [3] and rolling behaviors [4] have been reported. Couder and coworkers [5] evidenced that close but below 
the Faraday instability threshold, the bouncing of the drop becomes subharmonic, and Faraday waves are generated 
at each impact. The resonant interaction between the droplet and the waves emitted leads to a self-propulsion 
mechanism and the droplet becomes a walker. Spectacular phenomena are therefore obtained from this peculiar 
wave-droplet association. For example, Perrard et a/., using oil droplets with a ferrofluid core in a magnetic harmonic 
potential, obtained trajectories with quantified radii and angular momenta [6]. Using a circular coral, Harris et al. 
showed that the long term statistic of the walker dynamics looks like the PDF of an electron in a corral [ 7 ]. Each 
experimental works cited above considers 2d systems. Nevertheless, only a few investigations are based on droplets 
and Faraday waves evolving in a Id system [8]. As circular orbits can be considered as Id, one recent work focuses 
on strings of droplets propelled in annular cavities [9]. Confining droplet trajectories along Id systems is therefore of 
interest from a theoretical perspective and would provide a new experimental way to manipulate such droplets. In the 
present paper, we propose to study the behavior and dynamics of a single walking droplet within narrow submerged 
rectangular cavities, as sketched in Fig. [^a) and pictured in[^b). We will see that our objective to confine such 
droplet into Id motions will be reached and that the walker dynamics within the channels shows interesting analogies 
with electromagnetic waveguides. 


EXPERIMENTAL SETUP 



FIG. 1. (Color online) (a) Sketch of a rectangular cavity of width D = 2Xf- The liquid depth is adjusted in order to have 
h — A mm in the channel and ho = 1 mm elsewhere. The droplet can only walk in the channel as explained in the main text, 
(b) Picture of a walker above a channel of width D = 2Xf. Dashed lines indicate the channel boundaries. Faraday waves 
excited by the droplets are seen to be strongly damped outside the cavity. 


A container of section 180 x 180 mm^ is filled with silicone oil having a kinematic viscosity u = 20 cSt, a density 
p = 949 kg/m^ and a surface tension a = 20.6 mN/m. We use two types of containers: a square cell used for 2d 
comparisons and a cell with channels of 120 mm length but of different widths carved into it. We use a dispenser to 
generate on-demand droplets m, with a selected diameter 2R = 0.8 mm and a precision of 1 % . The tank is vibrated 
sinusoidally and vertically, using an electromagnetic shaker with a tunable amplitude A and a fixed forcing frequency 
/ = 70 Hz. An accelerometer is fixed on the vibrating plate, and delivers a tension proportional to the acceleration in 
order to determine the dimensionless forcing acceleration, as defined by F = A'K^Ap/g^ where g is the acceleration of 
the gravity. The acceleration F defines the so-called “Memory time” m tm = tfT i7/(ri7 — F) which measures the 
time of persistence of the waves emitted by the walkers. To ensure purely vertical oscillations of the container and to 
avoid parasite effects on the walking dynamics of droplets, we use air cushions carriage that surround the vibrating 
axis, and guide its vertical motion similarly to Ref. m- 
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WALKING REGIME 


The droplet is bouncing above a threshold T b [DlSllll]. By increasing the acceleration above a second threshold 
Tw > r^, the droplet becomes a walker that emits localized Faraday waves of wavelength Xp although the acceleration 
is below the Faraday threshold F p miis]. Walkers are only observed in the interval [rw,ri?]. Experiments made 
by Eddi et a/. [16] and Carmigniani et al HZ] show that modifying the depth of the fluid in the walking droplet 
experiment changes the value of the Earaday instability threshold F^? as well as the walking threshold Fw for a given 
frequency. In Fig. [^a), we plot the evolution of both thresholds Tp and F^v as a function of the liquid depth h 
for / = 70 Hz, delimiting three regimes: Bouncing, Walking and Faraday. It should be noted that this plot can be 
drastically modified by changing the forcing frequency /. Actually, we choose to work at / = 70 Hz to obtain the 
largest possible interval Tp — Tw for large liquid depths. This plot is quite instructive because below a critical value 
noted he ~ 1.5 mm, no walker is observed. However, one notices the apparition of the Faraday instability even for 
a small liquid depth but for large values of F. Finally, the cases with h > 3 mm correspond to a deep liquid region: 
F p and Tw are nearly independent of h. In our experiments, we fixed the forcing acceleration such that F = 0.95r p 
in the deep liquid regime, corresponding to a horizontal dashed line in Fig. [^a). There, the number of impacts still 
emitting waves is given by the memory parameter Me = ruj^p = 20. 




FIG. 2. (Color online) (a) For / = 70 Hz, thresholds of the Faraday instability Tp (triangle down) and the walking regime 
Fw (triangle up) as a function of the liquid height h. Three regimes are distinguished: Bouncing, Walking and Faraday. The 
horizontal dashed line corresponds to the experimental conditions (Me = 20) and circles indicates the two depths used in this 
work, (b) Thresholds of the Faraday instability (triangle down) and the walking regime (triangle up) as a function of the 
non-dimensionalized channel width D/Xp, herein the liquid height is fixed at /i = 4 mm. 


From the above results, we are able to design cavities ensuring that the droplet remains inside. Several channels 
were created, each one consists of a deep liquid area h = 4 mm surrounded by a shallow liquid area ho = 1 mm, as 
illustrated in Fig. [^a). These depths are measured with an uncertainty of ±0.03 mm. For a liquid depth h = 4 mm, 
we measured that the Faraday wavelength is Xp = 5.3 ±0.1 mm. Different channels have been created with a fixed 
length L = 120 mm and different widths D. Widths from D/Xp = 1/2 up to 6 have been used. In all cases, the 
droplet is confined in the cavity. In Fig. [^b), one can observe a walker evolving in a channel. The waves emitted by 
the successive bounces are seen to be strongly damped outside the rectangular cavity. Indeed, the Faraday threshold 
is far above the experimentally imposed forcing. In Fig. [^b), we plot the dependence of both thresholds as a function 
of the channel width. We note that, for large channels D/Xp > 3, Fp and Fw do not vary with the channel width, 
and their values correspond to those obtained for a square cavity in the deep water region, as shown in Fig. ia). 
In the experiments performed in smaller cavities: D/Xp < 3^ to ensure a memory parameter of Me = 20, one has to 
consequently adjust the forcing acceleration according to the evolution of the Faraday threshold within its repective 
channel. In our experimental protocol, F is tuned so that Me is kept constant. 
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WALKING IN CAVITIES 


Let us report various trajectories for droplets in different channels. For D/Xp < 1, the channel is able to pin 
the droplet which remains immobile, only bouncing on the surface. Walking droplets are only found for D/Xp > 1. 
Typical trajectories are drawn in Fig. for D/Xp ^ {1.5, 2,4, 5.5}. In narrow channels, the walker is following a 
rectilinear back-and-forth motion at a constant speed. When the droplet reaches an extremity, it is reflected. For wider 
channels, the situation is completely different since the droplet jiggles. Tortuous trajectories with several reflections 
on the sides of the cavity are seen. By increasing the width, the x component of the droplet velocity becomes more 
and more important. It should be equivalent to the y component for an isotropic and large system. 



FIG. 3. (Color online) (Top row) Typical trajectories of a walker in four channels: D/Xp ~ 1.5 (a), D/Xp — 2 (b), D/Xp ^4 
(c), D/Xp ~ 5.5 (d). In the first two narrow cavities (a) and (b), the walker follows a linear path. The droplet has a back and 
forth motion with a constant speed along the y-axis. In the wider cases (c) and (d), the walker is wobbling and oscillations 
are seen in the transverse direction, i.e. in the x direction. (Bottom row) Pictures of the Faraday pattern obtained in the 
corresponding upper cavities when the acceleration F reaches F p. One can observe evanescent waves outside the channels. The 
Faraday periodic pattern is only observed along the y axis. The pattern corresponds to a bump along the x axis, where the 
Faraday wavelength slightly differs from Xy. A secondary light source placed close to the liquid bath evidences a periodic 
substructure along the x axis, emphasized by a yellow dashed vertical lines. 


Since the successive impacts of the bouncing droplet on the surface are exciting Faraday waves, we explored the 
Faraday pattern into the different rectangular cavities. The forcing acceleration was cranked up until the emergence 
of Faraday waves. Pictures at the onset of the Faraday regime are shown in Fig. Instead of square patterns found 
in large tanks, we observed line patterns oriented along the transverse direction of the narrow cavities. The line 
spacing along the longitudinal direction corresponds roughly to Xp. Those patterns look similar to those obtained 
by Pucci and coworkers in the case of worm-like floating droplets [HHS]. We conclude at this point that a single 
mode of wave propagation could exist in narrow channels since only transversal structures are observed for F close 
to Tp. Outside the channel, in the transverse direction, one can notice evanescent Faraday waves, which are strongly 
damped, according to our results of Fig. [^a). Strictly considering line wave patterns, it is hard to understand why 
droplet starts to wobble for wide channels. A secondary light source has been placed close to the experiment in order 
to visualize the substructure of nodal lines which are evidenced with yellow dashed lines in Fig. One can notice 
that the second light has to be tilted enough to obtain a proper visualization of the nodal lines inside the channels. 
They are undulated with a wavelength close to along the transverse direction. Please note that A^, is slightly 
different from Xy We note a remakable feature: the number ni of nodal lines is found to be odd in those channels. 
This implies that an integer (and not half-integer) number of spatial wave oscillation periods is always realized in 
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our experiments. This peculiar property is exploited and discussed in more detail in the subsequent section. The 
number of nodal lines will affect the transverse speed as shown below. In order to quantify the walker dynamics 
in channels, we consider the average speed along the longitudinal and transversal directions: {\vy\) and {\vx\), as a 
function of the channel dimensionless width D/\p. The results are shown in Fig. |^a) and ib). One observes that, 
along the ^-direction, the velocity is equal to zero for channels with D/Xp < 1. As the width increases, the velocity 
suddenly grows to reach, for D/A > 3, a constant value, namely the velocity of a walker in a 2d cell without boundary 
effects. The behavior is different along the x direction. (|ua,|) has an almost constant value forl<T)/Ai?<2 while 
it steadily increases with D for D/Xp > 2. One expects to obtain {\vy\) = {\vy\) for infinitely large channels since 
one recovers, in this case, the 2D walker dynamics. The evolution of {\vy\) with D is described later in this article 
thanks to a electromagnetic wave-guide analogy. The evolution of (|ua,|) seems to be linked to the Faraday instability 
as discussed before. Indeed, when only one nodal line is observed (i.e. n/ = 1, which corresponds to one wavelength 
in the transversal direction), (|ua,|) is constant and only changes when two new nodal lines appears (i.e. ni = 3, which 
corresponds to the emergence of a second wavelength along the transversal direction) for D/Xp = 3. Note that the 
subsequent growth of {\vx\) is slightly less pronounced for2<D/Ai7<3 than for D/Xp > 3. We attribute this to 
the fact that in this particular interval of channel widths a single nodal line along the center of the channel is still 
most often encountered, while the theorical model developed in the next section would in principle allow also for the 
excitation of the subsequent mode exhibiting three nodal lines. This type of restriction is expected to become less 
effective for wider channels with D/Xp > 3. 

Since the main application of narrow cavities is the path control of walkers, the relevant question concerns the 
choice of the optimal width D for limiting the speed fluctuations along the transverse axis. Figure presents in a 
semi-log plot the ratio of the average kinetic energies measured along the x and y directions, i.e. the ratio {Vy)/{vl)^ 
as a function of D/Xp. As expected, a peak is observed around 1.5 < D/Xp < 2.25. It should be noticed that the 
ratio reaches values above 100. One can consider these range as optimal widths to perform Id experiments with 
walkers, in our conditions. 



DIKp 

FIG. 4. (a) Evolution of the average longitudinal speed (lu^/l) and (b) average transversal speed (|ucc|) of a walker as a function of 
the dimensionless parameter D/Xp. The black dashed curve is a fit using Eq. §. Different colors denote different propagation 
modes noted m, defined in section Discussion. Colors correspond to those of the trajectories shown in Eig. 
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FIG. 5. (Color online) Ratio of the kinetic energies in a semi log-scale as a function of the dimensionless width D/Xp. The 
ratio reaches a maximum at approximately 1.5 < D/Xp < 2.5. For very large D/Xp values, we expect that the ratio tends to 
1 due to the equiprobability of the droplet motion in both directions. 
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DISCUSSION 

The above observations can be qualitatively captured by means of an approximate theory of the Faraday waves 
that are emitted by the bouncing droplets. In a confinement-free configuration, a single bounce of a droplet at the 
position ro = (xo,^o) gives rise to a periodically oscillating surface wave profile of the form p!3| [2T] 

C(r, t) = AJo{kF\Y - rol) cos(wi/2)e“*/'^'^ (1) 

which decays with a rate 1/tm that is inversely proportional to Me. We neglect here the occurrence of a spatial 
attenuation oc of the surface waves, which is encountered in droplet experiments [21]. As we argue in more 

detail in Ref. [22], the Bessel function of the first kind Jo arising in Eq. 0 can be understood as representing the 
imaginary part of the retarded Green function G(r, that is associated with the two-dimensional Helmholtz 

operator, that satisfies the equation 

^ ^ ^^ 

We then have 

Jo(fcF|r - rol) = -4Iin[G(r, Tq, kp)] 

= ^ j d'^k5{k^ - 

which, when being inserted into Eq. 0, indicates that only those waves oc exp['ik • (r — ro)] persist whose wave vectors 
k satisfy the Earaday wave condition |k| = kp = ^'KjXp. The Earaday wave profile expressed in Eq. 0 will be 
modified in the presence of sub-surface channels. We suggest that this modification can be approximately accounted 
for by introducing additional boundary conditions to the Helmholtz equation, see Eq. 0 , that the Green function 
has to satisfy. Inspecting the experimentally obtained surface wave profiles that are depicted in the lower panels of 
Eig. 3, we infer that the Earaday waves emitted by a droplet in a channel approximately satisfy Dirichlet boundary 
conditions at the lateral borders of the channel, yielding the condition G(r, ro,/^^?) = 0 at r = [x^y) with x = 0 and 
X = D ioi the Green function of the Helmholtz equation. We define the horizontal coordinate system such that the 
channel lies within {) < x < D. They furthermore appear to exhibit a nodal line at the channel’s centre at x = D/2 
as can be seen again in Eig. 3, which is presumably required to allow for a stable droplet motion within the channel. 
We can then expand the Green function within the lateral channel eigenmodes 

Xm{x) = ^ sin (4) 
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satisfying Dirichlet boundary conditions and {—d‘^/dx^)xm{x) = k‘^Xm{x) with km = ‘Im'KjD^ and exhibiting a 
node at X = 1^/2. This then yields a superposition of longitudinal waves of the form Xm{.x) eyi^{ik!^^\y — ^o|) with 
= \/k‘^F — k‘^ for 0 < m < D/Xp^ complemented by another superposition of evanescent waves of the form 
Xm{x) — yo\) with = \/k‘^ — k‘jp for m > D/Xp the influence of which can be safely neglected. As 

a consequence, surface waves cannot be sustainably excited within too narrow channels whose widths D < Xp are 
below the Faraday wavelength. This in turn inhibits the walking of droplets within such narrow channels, which is 
in excellent agreement with our experimental findings. It is now straightforward to infer that the walking speed Vy 
of the droplet within the channel is proportional to the effective wave number that characterizes the longitudinal 
Faraday wave pattern along the channel, as this particular wave number determines the slope of the surface wave 
profile through which the droplet is horizontally accelerated upon impact. In the particular case of channels whose 
widths satisfy Xp < D < 2Xp, this effective wave number is evidently given by ky^^ = as only the lowest 

channel mode xi can be excited in that case. This yields the longitudinal walking speed 

= (^) ( 5 ) 

where Vy^^ is a characteristic speed scale of the droplet without constraints, i.e. the speed of a drop in 2d. 

Equation ^ is fitted on the data of Fig. |^a), yielding as sole free parameter Vy^^ = 11.9mm/s, which roughly 
corresponds to the speed for droplets in 2d. The fit is in good agreement with our experiments, even for D > 2Xp. 
This is attributed to the fact that longitudinal speed of the droplet is predominatly governed by the largest variation 
of the Faraday wave along the channel, and the latter would be provided by the lowest transverse mode xi. The 
transverse speed on the other hand, should be determined by the highest transverse mode that enters into the 
superposition of modes constituting the Faraday wave profile, since this particular mode provides the largest wave 
profile variations across the channel. We therefore should generally expect \vx\ oc m where m = rrimax is the largest 
integer that satisfies m < D/Xp. This expectation is roughly satisfied as we can see in Fig. 4(b). There are, however, 
significant deviations due to the fact that in practice this highest transverse mode Xmmax hardly ever participates 
in the experimentally observed Faraday wave profile, as we saw in Fig. 3. A more elaborate theoretical framework, 
taking into account the decay of Faraday waves outside the channel in a more quantitative manner, is certainly needed 
in order refine the model under consideration and thereby obtain a more detailed understanding of the behaviour of 
walking droplets in channels. The properties of the walkers dynamics reminds the propagation of light in rectangular 
electromagnetic wave-guides. In particular, one may notice the “Ray-Optics approach” used to understand the wave 
propagation in a guide, where the rays of light bounces on each wall. Let us compare both systems: wave-guides 
are usually made of two materials with refraction index ni and n 2 . The light propagates in the medium of index 
rii <722. Therefore, the light is focussed in the guide thanks to total reflection. In the walking droplet experiments, 
the two “materials” are the regions of different fluid depths which alter the propagation of the Faraday waves, as seen 
in Fig. The walker is constrained within the channel thanks to the differences between the Faraday thresholds. 
More interestingly, the analogy can be pushed further by comparing the velocity of the walker to the group velocity 
of the light propagating within an electromagnetic guide. One can compare here the velocity of the walker to the 
group velocity of the light propagating within an electromagnetic guide. Indeed for TEl wave guides, the light group 
velocity is given by Eq. (§. Nevertheless, some differences occurs. Indeed, the Earaday instability auto-adapts to 
the channel width where, in the electromagnetic case, the wavelength is an experimental fixed parameter. Therefore, 
the comparison breaks down for large channels which should be compared to TE^ wave guides. Eurthermore, this 
waveguide analogy does not give any explanation for the evolution of {\vx\) with the width of the channel. Einally, the 
wave emitted by the droplet are standing waves while electromagnetic wave-guides consider only propagating waves. 
Analogies should therefore be taken with caution. 


SUMMARY 

In summary, we evidenced that it is possible to confine and transport walking droplets along Id channels. While 
the longitudinal motion is dominated by a single mode, a fine structure composed of m modes is observed in the 
transverse direction. We discussed the similarities and differences with waveguided systems. We have shown, thanks 
to an energetic study, that the optimal width maximizing the longitudinal speed is around 2Xp such that applications 
can be designed. In particular, annular cavities have been investigated where droplets share a common coherent wave 
that propels the group at a velocity faster than a single drop [9]. In addition, upon shedding light on the following 






question: how straight the walker’s motion can be in a Id channel, this study completes the works based on walkers 
in confined geometries. Among them, one can note the numerical study of Gilet [20] and the experimental work of 
Harris et al [7] where the authors focused on the statistics of walkers in a corral cavity and cavity modes. We plan 
also to study branched structures in order to manipulate walkers and organize collisions or self-interferences. Such 
narrow cavities are therefore the beginning of many original investigations. 
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